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Tracking the optimal speed profile in electric train operation has been 
proposed as an efficient and feasible solution for not only reducing energy 
consumption, but also no at costs to upgrading the existing railway systems. 
This paper focuses on finding the optimal speed profile based on 
Pontryagin’s maximum principle (PMP) while ensuring the fixed running 
time, and comparing energy saving levels in the cases of applying or not 
applying PMP. The way to determine the fixed running time also differs 
from works published is to calculate the total trip time equal to scheduled 
timetable exactly. Calculating accelerating time ta, coasting time tc, braking 
time tp via values of maximum speed v», braking speed v» of optimal speed 
profile. The other hands, va and vp are determined by solving nonlinear 
equations with constraint condition: the running time equal to the demand 
time. Simulation results with data collected from electrified trains of Cat 
Linh-Ha Dong metro line, Vietnam show that energy reduction for the entire 
route when PMP utilization is up to 8.7% and running time complied with 
scheduled timetables. 
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1. INTRODUCTION 

High levels of air pollution and congestion are knotting issues related to transport in big cities. 
Under the circumstance, urban railway is considered a great potential to be efficient, reliable, and 
environmentally friendly means of transportation. However, its major disadvantage lies in consuming a large 
amount of electric energy in operation while the energy prices are steadily rising [1], [2]. Therefore, energy 
saving in train operation has become a hot topic to draw more attention from researchers worldwide. There are 
some energy efficiency strategies having been proposed. The first approach is related to enhance regenerative 
braking energy recuperation which can reduce the net energy consumption approximately 30% [1], [3], [4] by 
onboard/stationary energy storage systems [5]-[8], equipping the traction substations with reversible converters 
or active rectifiers so as to pump back the regenerative braking energy into utility source, as a result, all 
regenerated energy can be recovered [9], [10], and optimizing scheduled timetables makes the regenerated 
energy among trains interchange easier [11], [12]; the second one is energy-efficient traction systems, 
determining optimal position for substations by reducing energy losses in the power supply network, in on-board 
traction equipment [4]; the third one is energy - efficient driving by applying optimal theory is to minimize 
operating energy [13]-[18]. Among above approaches, the energy-efficient driving solution does not need to 
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invest in equipment or infrastructure of existing metro lines, so this one is most suitable for metro lines in 
Vietnam to have just been installed. The researchers of South Australia University — Albrecht et al [19] and 
Howlett et al. [20] outlined an energy-efficient driving strategy for a train journey on a track with uphill and 
downhill gradients by designing control laws to calculate location of optimal switching points, and then 
determining the optimal speed profile, Baranov et al. [21] proposed a solution to minimize energy 
consumption and consider fixed trip time by supplementing Lagrange multi-plier 2 in objective function, 
but calculating A to find the actual time equal to demand time is not easy. In addition, some researchers have 
also applied intelligent train operation algorithms, such as: using fuzzy logic, neural networks, genetic 
algorithms, and particle swarm. To improve its performance; Lu et al. [22] applied three optimization 
algorithms (colony optimization algorithm, dynamic programming, genetic algorithm) to search for the 
optimal speed trajectory, Acikba and Söylemez [23] utilized artificial neural networks and genetic algorithms 
to determine the optimal coasting point. However, these studies also did not mention fixed running time. 
For the continuous control problem, the Pontryagin’s maximum principle (PMP) is used to find necessary 
conditions for an optimal strategy. It is shown that these conditions yield key equations that determine the 
optimal switching points. Therefore, in this paper, PMP has been presented to determine optimal speed 
profile for Cat Linh-Ha Dong metro line in Vietnam, ensuring fixed - running time as well. Simulation results 
with different scenarios are showed the effectiveness of optimal control method - PMP in saving energy of 
train operation without any changes about infrastructure, facilities. 


2. TRAIN MODELING 
The urban electric train often operates in three motion regimes: accelerating, coasting, braking as in 
Figure |. To build-up the kinematic equation of the train, it is necessary to analyze the forces acting on the 
train, surveying the longitudinal profile of the train. Thereby, its kinematic equation can be represented by 
the following continuous - space model [24], [25]. 
dx 


dt as (1) 
mv Fy (v) — For) — Wo(v) — Forad (x) 


Where v,t,x,m represent respectively train speed (m/s), operation time (s), train position (m), 
full load translating mass of train (tone) and Fir, Fpr, Wo, Fgraa are traction, electrical braking, resistance, 
gradient resistance forces applied on the train. Depending on the structure of the trains, different gauges will 
have different impact forces and influence coefficients. Forces used in this paper are calculated from Cat 
Linh-Ha Dong metro line. 
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Figure 1. A typical speed profile with three motion phases for short inter-stations 


Base on curves of traction force F;,, braking force F,, given by manufacturers [26], using the 
identification method: least square method to find equivalent polynomials. The maximum traction and 
maximum braking forces corresponding to the speed v are: 
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é -í 13.2 (0 < v < 32) 5 
tr — (-2.5 - 1075y? + 0.007 - v? — 0.66v + 28.35 (32 < v < 80) @) 
14.7(0 < v < 65) 
Fy, = $ —0.254v + 31.21(65 < v < 75) (3) 
—0.2027v + 27.36(75 < v < 80) 


The basic resistancewgcan be calculated by using Davis formula [27]. 


Wo = 2 = at but cv? (4) 


Where a, b,c are coefficients of train’s resistance. The gradient force F,,,q caused by slope of road. 
Foraa = Mgsina (5) 


In which g, @ are the gravity acceleration and the rail track slope respective. 


3. PONTRYAGIN’S MAXIMUM PRINCIPLE APPLICATION FOR OPTIMAL SPEED 
PROFILE DETERMINATION AND FIXED RUNNING TIME 
There are several optimal control theories applied to find the optimal speed profile. In this section, 
Pontryagin’s maximal principle is applied to detect the optimal switching points, and then determining the 
optimal speed profile, ensuring fixed trip time. The main objective is to evaluate the energy savings when 
trains track the optimal speed profile compared with the original speed profile. 


3.1. Speed profile optimality analysis based on PMP 

Depending on the distance between stations, trains operate in four phases: accelerating, holding, 
coasting, braking. However, for urban metro lines, due to the short distance between stations, they only 
operate in 3 modes: accelerating, coasting, braking. From train’s state equations with boundary conditions, 
and the objective function J, using PMP seeks the optimal switching points of these modes; thereby 
determining the speed, distance, and trip time in each phase. 


3.1.1. Problem formulation 
The motion of a train along a track is rewritten from the state equations (1). 
dt 1 


a 2 (6) 
ara Utr ftr V) — Uprfor(V) — Wo) — fgraa x) 


Where Utr, Upr are defined traction and braking control variables of train, both of which are restrained by: 
Utr E [0,1]; upr € [0,1];Um» € [0,1]; fers fors fgraa are forces per unit mass; traction force applied at the 
wheels, braking force, mechanical force, gradient force acting on the train. 

Where: 


_ F 


Fp Fgrad 
fr = E; for = F; foraa = L (7) 


Therefore, boundary conditions are given by: 


v(0) = 0, v(X) = 0,t(0) =0 ; 
ee e eae (8) 


Where V(x) is the maximum allowable speed, X is the terminal of the train operation; v(0), v(X) are the 
speed at the beginning, at the end of the route; T is duration of the trip is also given by the timetable. 
From state equations (6), with boundary conditions (7), (8) how to find objective function in order to 
minimize the train’s operation energy consumption as the train runs from location x = 0 to location x = X in 
time T by controlling the traction force, while ignoring electric braking force since regenerative braking 
energy is not recovered. The objective function is written as: 


J = Sy ter fer )dx > min O) 
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3.1.2. Solution 
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By PMP finding optimal solutions of an objective function is equivalent to maximizing its 


Hamiltonian equation. Based on (6), (8), a Hamilton function is formed as: 


Uerftr(V)—Uprfpr()—Wo graa N 


v 


H = -uy fr) + Pr= + Po ( 


Where p1, p2 are adjoint variables. 
Adjoint variable differential equations are reformed: 














dpi _ OH _ 
dx ot 
dp2 ðH Oftr E 
T E By T ta gg + Paga + pa lrfe) = unr fir(¥) = wolo) = faraa 
_ P2 | Oftr _ Ofpr _ 9Wo 
vl av br av dv 
Define 
— P2 Seabee dp2 _ d(p-v) __. dv dp 
p=" sop v = p2. Therefore ae Pa ae 
=> @ = 1 (82 _ 2%) 
dx dx dx. 
Given: 


dv _ Utrftr©)—Uprf prV)—Wo)-f graa®) 


dx v 


Therefore, Hamiltonian function is rewritten. 
p 
= (p — D)utr ftr — Pupr for — P(Wo + foraa) + me 


Hamiltonian function is maximized by the following values of u, and upr: 


Utr = 1 if p>1 Upr = 0 if 0<p<1 
Ur E (0,1) if p= 1and4upr E (0,1) if p=0 
Utr = 0 if p<1 Upr = 1 if p<0 


From the above analysis, five optimal control laws are designed: 
- Full power (FP):u;, = 1,up, = 0 when p > 1 
- Partial power (PP):u;, € [0,1], upr = 0 when p = 1 
- Coasting (C): uz, = 0, upr = O0when0<p<1 
- Full braking (FB): uz, = 0, upr = 1 when p < 0 
Partial braking (PB):utr = 0, upr € [0,1] when p = 0 
Subatients (13), (12) in (14), finding the differential equation for p(x) 


dp _ (1-p) 
dx v 





Utr fir) +- © Uprfor(V) +? zwo w) -35 


From (11), easily, p,is chosen by 0. 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


Full power mode: p > 1, up, = 0, Utr = 1, finding accelerating time tą, accelerating distance xg using (18). 


dp _ (1-P) pı P 
a a; feU) + wow) 


(19) 


From (6) finding the differential equation to determine Xa, ta, with initial conditions: x(0) = 0, t(0) = 0. 








dx _ v 

dv Uerf trV)—Wo(Y)-f grad(*) 
dt 1 

dv Uerf trV)—wo(Y)-fgraa(*) 
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Coasting speed vņis calculated as following [20]. 


_ Wn) 
~ o'n) CD 


Vp 
Where vp - hold speed: 
~(v) = v-w(v) = v(a + bv + cv’), W(v) = v? - wo(v) = v? (b + 2cv) (22) 


From (6) finding the differential equation to determine Xe, te, with t(v = Vp) = ta, x(V = Vn) = Xa 








dx _ v 

dv —wo(v)-f grad(*) 

ee (23) 
dv wo(v)+f graa(*) 


Full braking mode: p < 0, up, = 1, Utr = 0, finding braking time tp, braking distance x, using (14). 


d 


E= E for) + Eww) (24) 


From (11) finding the differential equation with t(v = vp) = tp, x(V = vp) = Xp 





dx _ v 
dv —Uprf br(Y)—Wo(Y)-f graa(*) 
dt 1 (25) 


dv Uprfpr(Y)—Wo(v)-F graa(X) 





3.2. Fixed running time 

To meet requirements of train’s running time and stopping time at stations as the scheduled 
timetable when trains track the optimal speed profile. In this section, we will calculate t, - the time of the 
accelerating phase, t, - the time of coasting phase, t, - the time of braking phase, and the total running time 
the of three phases must abide by the scheduled timetable. Calculating ta, tp, to is conducted by Maple 
software. 


3.2.1. Accelerating phase 
Equation motion of the train in optimal traction mode: 


dv 
PR fir(V) — Wo(v) (26) 
Using the variable dissociation method, the running time in accelerating phase is expressed as: 


ae dt > f*—_" __= f*dts t — st (27) 


Fir(®)—wo(v) 0 Fer(v)—wo(v) 0 Fr@—w%  *¢ 


Where: fis calculated as (7). 
From (27) the acceleration time may be employed: 


32 Vh 
dv dv 
o fer@™-wolv) Jy, Fer)-wo(v) 


32 
v 
_ dv +Í k dv 
=< 13.2 2 —5 3 2 2 
— —2.5:10 +0.007-v“—0.660+28.35)—(a+bv+ 
h 247000 (at+bv+cv2) 32 ( v v v )-(a+bv+cv*) 








32 (28) 


1 
¥1014ac—25-1012b2-5.34-109¢C 








-1013 
arctan (250 -107? 2:10°"(2cv+b) ) 


V¥1024ac—25-1012b2-5.34-109¢ 
vh 





0 





15.25:1015 ac—3.81:1015b2—31.89-1012 a—936.16-1012 b— 


-izssaotareran( 
33.72-:1015c+13.04-1012 


sag ea pe ee e) 





32 
¥15.25-1015ac—3.81-1015b2-31.89-1012a—936.16-1012b—33.72-1015c+13.04-1012 
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Using the Maple software tool, the acceleration time as a function of velocity. Acceleration distance is 
calculated as: 


dx vdv x v vdv v vdv 

Z = v > dx = vdt = ————__> f" dx = |? ——_ > x, = [M 29 

dt ftr@)—wo(v) So So ftr(v)-wo(v) S So fer&)—wolv) @) 
3.2.2. Coasting phase 

Motion equation of the train in optimal coasting mode: 

dv 

= = —wo(v) (30) 


Using the variable dissociation method, the running time in coasting phase obtains: 











Vb Vh 
dv dv tatt dv 2 2cv+b 
= dt > = |A “dtot, = tg = arctan ( ) 31 
-wo (v) vh —Wo(v) Sea S vp a+bv+cv? /4ac—b2 V4ac—b? GD 


In which braking velocity v, is given as: 








Vp = ZZN withp(v) =v: w (v) = v(a + bv + cv’), W(v) = v? -wi (v) = v? (b + 2cv) 
h 
_ Yh a Vh 
Vp = 1400h) ~~ atvp(b+cvp) (32) 
VAWo(h) vp(b+2cvp) 
Coasting distance is computed as: © Z v > dx = vdt = 
dt -wo (v) 
XatXe _ [Yb vav _ [Yh vdv 

L dx = Vh -wo (V) kes i at+bv+cv2 (33) 
3.2.3. Braking phase 

Motion equation of the train in optimal braking mode: 

dv 

2 = —wo(v) — fr(v) (34) 


Using the variable dissociation method, the running time in braking phase is given by: 


dv 0 dv tattc+t v dv 
—= dt> | — = [h gt > [P — Z =t 35 
—for)—-wol(r) Sin —forv)—wo(v) tatte So for(@)+wo(v) k (23) 


Fy is calculated as (7), from (35) the braking time can be written as: 


Vb 
dv 1 2cv+b 
t=] 1—2 = 2 arctan (1 (36) 
0 fbrmaxt(atbvt+cv?) 4ac—b*+4cf brmax  4ac—b*+4cf prmax 


vdv 


The braking distance is calculated: from © Z v > dx = vdt = — _ 
dt —fbr(v)—wo(v) 


XatXcetXp _ r° vdv _ cb vdv 
i +X, rs Sy - - HAET Í 
atXc b -f pr(v)-wo (v) for)two) 


(37) 


From the above calculations, total trip time is equal to tg + tp + te 


4. SIMULATION RESULTS 

The simulation is based on the data of Cat Linh-Ha Dong metro line, Vietnam. There are 12 stations, 
1 depot, 6 traction substations, and two-side power supply mode. In this paper, simulation results are 
performed for the first Cat Linh station to the 12th Yen Nghia station with 12.661 km in length [28]. 
Simulation parameters are shown in Table 1 and Table 2. Cat Linh - Ha Dong metro line is short, operation 
modes of electrified train are comprised of three regimes: accelerating — coasting > braking. 
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Table 1. Simulation parameters of train 

















Parameters of metro train Unit Value Parameters Value 
Train grand-up 2M2T a 1.19 - 107? 
Full load translating mass kg 246700 b 2.56- 107° 
Number of electrical traction units 08 c 1.54- 1074 
Max speed km/h 80 
Base speed km/h 40 
Dwell time s 30 
Max acceleration/braking rates m/s 0.94/1 





Table 2. David’s coefficients of train’s resistance 


Figure 2 shows the responses of the optimal speed curve over the whole route in red when applying 
Pontryagin’s maximum principle, and the speed curve without applying this principle in blue. The optimal 
speed profile has lower accelerating speeds, and longer coasting distances. Therefore, less power is 
consumed on the optimal speed profile. The responses in Figure 3 also indicate that the optimal speed profile 
and the original speed profile ensure the same running distance between stations and the fixed trip time 


according to the schedule. 
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Figure 3. A comparison of optimal time profile and original time profile 
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Table 3. Results of a comparison of energy consumption with / without energy optimal strategy 








F Maximum Practical Maximum speed Optimal 
; Trip ; 5 ; Energy 
Inter-station length Distance fee Speed in energy in optimal speed energy saving 
(m) (s) original speed consumption profile consumption (%) 
profile (km/h) (kWh) (km/h) (kWh) 
Cat Linh - La Thanh 931 88 54.5 8.31 51.92 7.50 9.75 
La Thanh - Thai Ha 902 78 60 10.20 57.75 9.40 7.84 
Thai Ha - Lang 1076 91 60 10.20 59.06 9.86 3.33 
Lang - Thuong Dinh 1248 103 64 11.73 61 10.60 9.63 
Thuong Dinh - Ring Road 3 1010 79 68 13.41 65 12.23 8.80 
Ring Road 3 - Phung Khoang 1480 104 75 16.75 713: 15.82 5.55 
Phung Khoang - Van Quan 1121 86 69 13.85 66.27 12.66 8.59 
Van Quan - Ha Dong 1324 97 73 15.74 69.7 14.17 9.97 
Ha Dong - La Khe 1110 84 70 14.30 66.27 13.18 7.83 
La Khe - Van Khe 1428 101 75 16.75 72.58 15.53 7.28 
Van Khe - Yen Nghia 1032 81 68 13.40 64.77 12.04 10.15 
Total energy consumption 144.64 132.99 8.7 





Regarding as track conditions, constraints, the speed from a station to another station is different. 
Results show in Table 3 the slowest speed at 54.5 km/h, the highest speed at 75 km/h, and being always 
smaller than limit speed 80 km/h while trip time complied with scheduled timetable. Table 3 also 
demonstrated that the lowest saving energy is 3.33% while the highest saving energy is to 10.15%. Therefore, 
saving energy of the whole route is 8.7%. 


5. CONCLUSION 

Simulation results by applying Pontryagin’s maximum principle to find the optimal speed profile 
with fixed trip time for 12 stations of metro line Cat Linh - Ha Dong, Vietnam with three operation phases: 
accelerating, coasting, braking showed saving energy percent is up to 8.7% comparison with the original 
speed profile while the running time is unchanged. This research has not only brought remarkably economic 
benefits in energy consumption reduction without any costs for infrastructure of existing metro lines, but also 
created a foundation for designing the optimal speed profile applied to automation train operation (ATO). 
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